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In the last few years, several approaches have been developed to compute the exact Casimir 
interaction energy between two nonplanar objects, all lead to the same functional form, which is 
called the TGTG formula. In this paper, we explore the TGTG formula from the perspective of 
mode summation approach. Both scalar fields and electromagnetic fields are considered. In this 
approach, one has to first solve the equation of motion to find a wave basis for each object. The 
two T's in the TGTG formula are T-matrices representing the Lippmann-Schwinger T-operators, 
one for each of the objects. Each T-matrix can be found by matching the boundary conditions 
■ imposed on the object, and it is independent of the other object. However, it depends on whether 

the object is interacting with an object outside it, or an object inside it. The two G's in the TGTG 
formula are the translation matrices, relating the wave basis of an object to the wave basis of the 
, other object. These translation matrices only depend on the wave basis chosen for each object, 

and they are independent of the boundary conditions on the objects. After discussing the general 
I theory, we apply the prescription to derive the explicit formulas for the Casimir energies for the 

sphere-sphere, sphere-plane, cylinder-cylinder and cylinder-plane interactions. First the T-matrices 
for a plane, a sphere and a cylinder are derived for the following cases: the object is imposed with 
Dirichlet, Neumann or general Robin boundary conditions; the object is semi-transparent; and the 
object is a magnetodielectric object immersed in a magnetodielectric media. Then the operator 
approach developed by Wittman [IEEE Trans. Antennas Propag. 36, 1078 (1988)] is used to derive 
the translation matrices. From these, the explicit TGTG formula for each of the scenarios can be 
written down. On the one hand, we have summarized all the TGTG formulas that have been derived 
CIh' so f ar f° r the sphere-sphere, cylinder-cylinder, sphere-plane and cylinder-plane configurations. On 

the other hand, we provide the TGTG formulas for some scenarios that have not been considered 
i-C ' before. 
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Casimir effect is a purely quantum effect that arises from the vacuum fluctuations of a quantum field in the presence 
of boundaries It has attracted a lot of attention due to its wide applications in different areas of physics such 
as quantum field theory, condensed matter physics, atomic and molecular physics, astrophysics, gravitation and 
cosmology, nanoscience and mathematical physics |2|. A number of books and reviews have been written on this 
subject [ItDl. 

Naively, the Casimir energy is a sum of the ground state energies of all the eigenmodes of the quantum field. 
Various approaches have been developed to compute the exact Casimir energy analytically, such as mode summation 
method, path integral approach, Green's function method, quantum statistical approach and scattering approach. 
Nevertheless, before the turn of the century, the exact computations were limited to some simple configurations, such 
as two parallel plates, a sphere, a cylinder, two concentric spheres, a wedge and a cone. A particularly large amount 
of works were devoted to two parallel plates. 

Since the work of Lamoreaux [Tij in 1997, Casimir force have been measured with high precision in various configu- 
rations especially the sphere- plane configuration [p], EH ■ This has called for precise computation of the Casimir force 
between a sphere and a plate in particular, and between any two objects in general. Several approximation schemes 
have been developed for this purpose, such as the semiclassical approach [la - [l7( . the optical path approximation 
[l8l - |20j and the multiple reflection approximation [2l|. However, each of these numerical methods have its limitations. 
Therefore, it is compelling to have an exact formula for the Casimir interaction between two objects, so that one can 
justify the accuracies of the numerical approximations and the experimental data. 
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In 22 26], Gies et al. derived a worldline representation of the Casimir interaction between two objects imposed 
with Dirichlet boundary conditions. They used the Monte-Carlo technique to develop a numerical scheme to compute 
the Casimir interactions of various configurations from the worldline representations. This is basically a geometrical 
approach, and its application has so far been limited to scalar field with Dirichlet boundary conditions. 

Another series of breaktrough appeared in 2006. In 27], Bulgac et al. computed the Casimir interaction energy 
between Dirichlet spheres or between a Dirichlet sphere and a Dirichlet plate using the multiple scattering approach, 
whose application in Casimir effect can be dated back to the work of Balian and Duplantier 28, 29] . In this approach, 
the Casimir energy is expressed as an integral over the density of states of the fluctuating field, and the Krein's formula 
is used to relate the density of states to the scattering matrix of the objects. A novel feature about this method is 
that it is free of ultraviolet divergences. 

In (30j , Emig et al. used path integral quantization with partial wave expansion to compute the Casimir interaction 
energy between a cylinder and a plane when both are imposed with Dirichlet, Neumann or perfectly conducting 
boundary conditions. This was elaborated and extended to two cylinders i n 1311. Their approach can also be put in 
the framework of multiple scattering theory, as is more obvious in the works [32l - l36| . In , a general scheme that uses 
multipole expansions and multiple scattering theory to compute the Casimir interaction energy between two dielectric 
objects was proposed. This scheme was applied to obtain the Casimir interaction energy between two dielectric spheres. 
It was further illustrated in [33j for scalar interactions and in [36j j for electromagnetic interactions. The Casimir 
interaction energies between two spheres imposed with Robin boundary conditions, between two perfectly conducting 
cylinders, between a magnetodielectric sphere and a magnetodielectric half-space, and between a magnetodielectric 
cylinder and a magnetodielectric half-space were derived. 



Back in 2006, Bordag 37] has also used the path integral and Green's function method to find the Casimir interaction 
energy in the cylinder-plane and sphere-plane configuration for a scalar field with Dirichlet or Neumann boundary 
conditions. Later he also obtained the Casimir interaction energy between a semitransparent cylinder and a plane 
[38j | . The path integral formalism used in these works can be dated back to the works [39|, H(| ■ 

At the same time, Dalvit et al. j4l], IHJ used an entirely different approach - the mode summation method, to 
compute the Casimir interaction energy of two eccentric cylinders with Dirichlet, Neumann or perfectly conducting 
boundary conditions. Later in (43j, the case where the cylinders are separated by a dielectric medium was considered. 

In [iil . |45| , Kenneth and Klich explored the multiple scattering approach from the point of view of the Lippmann- 
Schwinger T-operator. In [iij . they have obtained a general formula for the Casimir interaction energy prior to [53 
and showed that the Casimir interaction force between two objects related by reflection is always attractive. In [45| . 
they used their formalism to derive the Casimir interaction energy between two dielectric spheres. 

In 0, |47j , Milton and Wagner used the multiple scattering approach to compute the exact Casimir interaction 
between two semitransparent spheres and two semitransparent cylinders, and obtained their weak coupling limits. 

Despite the different perspectives, all the methods, except for the mode summation approach of Dalvit et al., are 
fundamentally equivalent. In all cases, the formula for the Casimir interaction energy can be cast in the form 



E Cas = — d£1Hn(I-TTiGi 2 T 2 
*n Jo 



0, 



which is referred to as the TGTG formula. Different approaches lead to different methods for computing the Ti , T 2 
and Gi 2 ,G 2 i matrices. In multiple scattering theory, Ti (T 2 ) is the transition matrix or Lippmann-Schwinger T- 
operator of object 1 (object 2), which is related to and can be computed from the scattering matrix of the object. Gi 2 
(G 2 i) is the translation matrix that relates the wave basis centered at object 1 (object 2) to the wave basis centered 
at object 2 (object 1). 

As discussed above, the exact Casimir interaction energy between two objects have been mostly computed by the 
multiple scattering approach. Since mode summation have been proved to be a very powerful method for computing 
the Casimir energy of a cylinder, a sphere, two concentric cylinders or two concentric spheres, it is timely to extend 
the mode summation approach to Casimir interactions between two or multiple objects. This is the task undertaken 
in this paper, with the hope that one can get some new insight about how to compute the T and G matrices. 

In the first part of this paper, we explain how the mode summation approach can be used to derive the TGTG 
formula for the Casimir interaction energy between two objects, for scalar fields as well as electromagnetic fields. We 
consider both the case where the two objects are outside each other, and the case where one object is inside the other. 
We also discuss in Appendix [A] how this approach can be generalized to more than two objects. Along the way, we 
obtain prescription for computing the T and G matrices. 

In the second part of this paper, we illustrate the mode summation approach by considering the Casimir interactions 
between two spheres, two cylinders, a sphere and a plane, and a cylinder and a plane. First we compute the T-matrices 
for a plane, a sphere and a cylinder. Three cases are considered: the object is imposed with Dirichlet, Neumann or 
Robin boundary conditions; the object is semitransparent; and the object is magnetodielectric. As a matter of fact, 
scattering theory has been extensively used in classical and quantum field theories. Therefore, some of the T-matrices 
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(up to some constants) have been well-known and it is quite impossible for us to provide complete references here. 
We apologize for not citing any earlier works. 

The more technical part is the translation matrices. For two spheres or two cylinders, the translation matrices have 
been rather well-known for both scalar and electromagnetic interactions. For the sphere-plane and cylinder-plane 
configurations, the translation matrices for electromagnetic (but not scalar) interactions have been written down in 
[361 ] implicitly without explicit derivation, as the multiplication of a change of basis matrix and a translation matrix 
for plane waves. To give a unified treatment, we generalize the operator approach in [48j to find the translation 
matrices for the scalar and electromagnetic Casimir interactions of cylinder-cylinder, sphere-plane and cylinder-plane 
configurations. From these, we can write down explicitly the TGTG formulas for the scalar and electromagnetic 
interactions of each of these configuration. We obtain some new results that have not been considered before, 
such as the Casimir interaction energy between a semitransparent object and an object imposed with Dirichlet, 
Neumann or Robin boundary conditions; between a semitransparent sphere and a semitransparent plate; and between 
a semitransparent cylinder and a semitransparent plate. 

The layout of this paper is as follows. In Section [HI we present the general theory of mode summation approach 
for the derivation of the Casimir interaction energy between any two objects, when both are outside each other, 
and when one is inside the other. The electromagnetic case is discussed in detail. In Sections IIII1 IIVI and [V] we 
compute the T-matrices for a plane, a sphere and a cylinder, under various boundary conditions. In Section IVI1 
we compute the translation matrices for the sphere-sphere configuration and obtain explicit formulas for the scalar 
and electromagnetic sphere-sphere interactions. In Sections IVII1 IVIII1 IIX1 we do the same for the sphere-plane, 
cylinder-cylinder, and cylinder-plane configurations. 



II. GENERAL THEORY 



In this section, we interpret the TGTG formula for the Casimir interaction between two objects from the point of 
view of mode summation approach. We consider both the scalar fields and the electromagnetic fields. 
For scalar fields, the equation of motion is 



1^1 

C2 W 



I, 2 



+ V(x) )<p = 



(1) 



We consider two types of boundary conditions. 

In the first case, the potential function V(yi) is zero and the field <p> satisfies the boundary condition 



up + v 



dp 

On 



= 



(2) 



on the boundary surface S of the object. Here n is a unit vector normal to the boundary S and pointing away from 
the object. When v = 0, we can take u = 1, and this is the Dirichlet boundary condition. When u = 0, we can take 
v = 1, and this is the Neumann boundary condition. In general, if u and v are both nonzero, we can take v = 1, and 
this is the general Robin boundary condition with parameter u. 

In the second scenario, V(x) = A5(x) is a Dirac delta potential function with support on the boundary of the 
object. We say that the boundary is semitransparent. Integrating the equation (fTJ) across the boundary of the object, 
we find that the scalar field satisfies the boundary conditions 



ip — tp =0 
s + s_ 



dip 
dn 



dip 
On 



Xip. 



(3) 



Here S+ and S— denote respectively the outside and inside of the boundary surface S. 
Let 



/oo 
dw<p(x, w)e 
-oo 



We find that <£>(x, k) satisfies the equation 



V 2 </?(x, k) = -k 2 ip(-K,k), 



(4) 
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where the wave number k and the frequency u satisfy the dispersion relation 



k =^-—- (5) 

Fixing a coordinate system, we can express the solutions of the equation (j4} as a linear combination of the regular 
wave ip r ° s which is regular at the origin of the coordinate system, and the outgoing wave </?° ut which decreases to zero 
rapidly as x —¥ oo and k is replaced by ik. Here a labels the solutions. 
For electromagnetic fields, the Maxwell's equations read as 

V-D = P/ , VxE+— =0, 

<9D 

V-B = 0, VxH--=J f . 

at 

For a magnetodielectric object, the charge density pf and the free current density J/ on the object are zero, and we 
assume the linear relations 

D(x,w) =e(w)E(x,a;), 

H(x )W )=i-B(x, W ), (6) 

where e and \i are respectively the electric permittivity and magnetic permeability of the object. Here for Z = 
D,E,H,B, 

/OO 
dojZ(x,uj)e' iuJt . 
-oo 

The boundary conditions are the continuities of eE n , Ey , B„, —By across the boundary of the object. These conditions 

are not independent. The continuity of Ey implies the continuity of B„, and the continuity of —By implies the 

" fx " 

continuity of eE n . 

As usual one can introduce a vector potential A that satisfies the gauge condition V • A = 0. In terms of A, 

BA 

E = - — , B = VxA. 

at 

Then the Maxwell's equations are equivalent to the wave equation 

d 2 A 

V x V x A + £M-^2" = ( 7 ) 

Let 

/OO 
rfwA(x,w)e- ,wt . 
-oo 

The wave equation can be written as an eigenvalue problem, 

V x V x A(x, k) = fc 2 A(x, k), (8) 
where the wave number k and the frequency lo satisfy the dispersion relation 

k = y/sjiuj. 

The solutions of the equation ((5J) can be written in different bases. In the cases of rectangular, cylindrical and 
spherical bases, we can choose a distinct direction and divide the solutions into transverse electric (TE) waves A™ 
and transverse magnetic (TM) waves A™ parametrized by some parameter a and satisfy 

iv x A™ = A™, iv x A™ = A™. (9) 
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Moreover, the waves can be divided into regular waves A™' reg , A™ ,rcg that are regular at the origin of the coordinate 
system and outgoing waves A™' out , A™> out that decrease to zero rapidly when x — > oo and k is replaced by ik. 

Consider two objects Oi and O2. Choose appropriate coordinate systems with coordinate origins at O and O' 
respectively for each of the objects. Let x and x' be respectively the position vectors of a point with respect to O 
and O'. If L is the position vector of O' with respect to O, then x = x' + L. In the following, we derive the TGTG 
formula of the Casimir interaction energy between these two objects using mode summation approach. We consider 
separately the case where the two objects arc outside each other and the case where one object is inside the other. 
We will focus on discussing the case of electromagnetic fields. The case of scalar fields is easier and can be obtained 
in the same way. 



A. Two objects are outside each other 

Inside the object Oi, we express A in the coordinate system centered at O: 



A = f°° (ko^2 (^? A a E ' reg ( x >w) + CfA™' reg (x,w)) e~ iwt . (10) 
Inside the object O2, we express A in the coordinate system centered at O': 



A = 



/oo . 
dw£ te E ' re6 ( x » + ^A™' reg (x>)) e" 
-OO a 



-ibjt 

(3 

In the region outside the two objects, A can be expressed in both coordinate system. Close to the object Oi, 



/oo 
duj^ «A™' rcg (x,w) + 6?A™' out ( x ,cj) + c?A™' rog (x,w) + d?A™- out (x, w)) e~ iwt . (11) 
-00 a 

Close to the object O2, 

a f°° j ( a TE.rog/ / \ . ,fl A TE,out/ / \ , 8 » TM,rcg / / \ . A TM,out / / \\ 

A= / duj2_^U^A fj S (X ,LJ) +&2 A /3 ' ( x ^)+ c 2 A p (x,^)+rf2 A ( X )W)J. 

J —oo a 



In fact, the regular waves close to Oi are propagated from the outgoing waves outside O2, and the regular waves close 
to O2 are propagated from the outgoing waves outside Oi. More precisely, they are related by translations: 

/A TE '- ( x>)\ _vKJ TE (-D TE (-L)\ (A™~(x, W ) 

A ™,out (Xja;)v ) I f/ |feTM (L) ^,TM (L) I I A |M, reg(x>) I • 



From these, we find that 



|)=E^(-l)(| 

= X>, Q (L) (| 



(12) 



tr n • n o o , 1 ,• , • -, u 4 rrTE.TE rr TE,TM r ,TM,TE rr TM,TM 

Here D Qj( g is the 2x2 translation matrix with components U a ^ , U a p , l/ Q a , t/ a ^ 

Now, for i = 1,2, the boundary conditions on the boundary of the object Oj will give rise to four homogeneous rela- 
tions between af,bf,cf, df and Af,C". Eliminating Af and Cf , we are left with two relations among af,bf,cf,df. 
From these we can solve for bf , df in terms of af , cf and write it in the matrix form 

%) = - t * (i) • (i3) 
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Up to a constant, Tj is the transition matrix or the Lippmann-Schwinger T-operator associated with object Cv It is 
closely related to the scattering matrix of the object Oj 0, HU, [3(| |44L|451|. 
From (TT2l and (TT3l). we have 



T?^U Q ,^(-L) 

P I a 2 



=T?^U a ,^(-L)T 



2 „/3 



or in other words, 



53 ( I a , a / - 53 U Qlj8 (-L)1[f U^, a / (L) 



df7 V° 



Here I a ,a' is the 2x2 identity matrix. For nontrivial eigenmodes, the infinite vector with a-component given by 

(bf d") T cannot be identically zero, which holds if and only if det (I — M) = 0, where I is the identity matrix and 
M = T1U12T2U21. Tj is an infinite diagonal matrix whose (a, a)-component is the 2x2 matrix Tf, U12 is an infinite 
matrix with (a, /3)-component given by U Qi/ g(— L), and U21 is an infinite matrix with (/3, a)-component given by 
i]p a (L,). Therefore, we find that the eigenfrequencies of the system are the nonnegative solutions of the equation 

det (I-M(w)) = 0. 

Hence, the Casimir interaction energy is given by 

uj>0 
det(I-M(w))=0 

Using the residue theorem, this can be computed by 

£cas = 53 ^Res w -^- In det (I - M(w)) 

u>>0 

fc 1 />ioo 1 

= / dujuj— In det (I -M(w)) ( 14 ) 

/ ^e^Trln(I-M(^)). 
Jo d£ 



h 

2T 



Integration by parts give 



where 



E Cas ~ I day In III -r£[/;Vi. (15) 

Z7T 







M(tO - T 1 (iOUi 2 (iOT 2 (^)U 2 i(iO- 

Eq. ([T5|) is called the TGTG formula for the Casimir interaction energy. 
In the case of a scalar field, we let 



l p= I d W 53^ cg (x,c)e— * (16) 
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inside object Oi, 



/CO 
da;^Af^ es (x',a;)e--' 
-oo £ 



inside object O2. Outside Oi and O2, 

/oo 
dw 53 (of^(x,w) + 6?< ut (x, W )) 
-OO 



,00 f x (17) 

^ —00 a 



Then the Casimir energy is given by (|15[) . where the matrix T, is obtained from the boundary conditions on object 
Oj. It is diagonal with (a, a)-component T" being a scalar relating bf with af: 



bf = -T?a?. 



For the translation matrices U12 and U21, their components are defined by the relations: 

^ ut (x>) =E c/ ^(- L )^ es ( x ' w )> 
v r(x, W )=^(/ fta (L)^(x». 

B. One object is inside the other 

Assume that object Oi lies inside object O2. Inside the object Oi, we express A in the coordinate system centered 
at O: 

/OO 
dujJ2 (^?A™' rcg (x, u) + C? A™- rog (x,w)) e~ ia,< 
-00 a 

Outside the object O2, we express A in the coordinate system centered at O': 

A f°° J ( oP A TE,out / / \ . j-,/9 A TM.out / / \\ -iut 

A = / du^iB'^Ap (■x,L))+D%Ap (x,w)Je 

J-00 p 

In the region between Oi and O2, 

A/oo 
j \ * I a \ TE,rcg / , \ , 1a < TE.out / \ , a <TM.reg/ , \ , ict \ TM,out/ \\ — iwt 
^2^( tt i A « (x,w)+6iA a (x,w)+ Cl A Q ^x.tdj+^A,, ■ (x,w)je 
-00 Q 

which can be re-expanded with respect to O' as 

A /"°° J V* / ^ A TE , re S/ ' \ 1 ifi a TE.out/ / x . /3 * TM,rcg / / \ . ,/3 A TM,out/ / \\ -iut 

A= / dw2^( a 2 A /3 ( x > w ) + ^2 A /3 (x,w) + c^A^ S (x,w) + ^A (3 (x,w)Je 

J-00 ^ 

using the translation formulas 

/ A TE -(x>)\ _ v /<^ te (-L) F™-™(-L)\ (Al^ M \ 

[4^^'^) ~^kf ™(-l) Ur< rcg (x,-)J 



A^out^^x ^-,- (L) wi„ (l) \ a;— (x»' 

Ar' out (x, w )J-^ U™(L) W™-™(L)J Ia^V,"), 
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In other words, 

,0 



/3 

As before, the boundary conditions on the boundary of Oi gives 



-Tf Pi . (18) 



The matrix T is the same as in the case when two objects are outside each other. On the other hand, by eliminating 
?2 and I?2, 



£?2 an d f 2 j the boundary conditions on the boundary of O2 give rise to a relation 



= -n 1 ) • (19) 



Therefore, 



6?\ _ 



1 w 

-Tfx;va l/S (-L) (i) 
T?x;va,/ S (-L)^ Qj 



=T^V a ,^(-L)T^W^(L) ( b i 

P a' ^ 1 

As in the case where two objects are outside each other, one then obtains the formula (fT5|) for the Casimir energy, 
where now the matrix M is given by 



TiVi 2 (-L)T 2 W 2 i(L). (20) 



In the case of a scalar field, we let 

<P 

inside object Oi, 

<P 

outside object 2 . In the region between Oi and 2 



j — 00 



du^B^i*-',^-^ (21) 





- = l A,]T «^(x,c) + &?<C(x, W )) e** 



d,LU 



^(a^" s (x', W ) + &^^(x', W )) e - 



(22) 



The Casimir energy is given by (fl5|). with M given by (|20|1. The matrices Ti and T 2 are obtained from the boundary 
conditions on objects Oi and 2 respectively. They are diagonal with 
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For the translation matrices V12 and W21, they are defined by the relations: 

^ es (x',a;)=^T/^(-L)^(x,u;), 

(X 

^r(x,^=^^, Q (L)^ ut (x', W ). 

Using the mode summation approach, we have shown that the Casimir interaction energy between two objects can 
be expressed as the TGTG formula (|15[) . Our approach here is somehow formal in contrast to the mode summation 
approach employed in (4ll . |42| for the Casimir interaction between two eccentric cylinders. The Casimir self energies 
of each of the objects do not appear in our approach, and we obtain immediately the Casimir interaction energy. 
Nevertheless, we would like to emphasize that our approach can also be made rigorous as in 

mm- 

In the above and in the following, the Casimir energy is understood as the Casimir energy at zero temperature. We 
will not elaborate on the finite temperature Casimir free energy as it can be easily derived from the zero temperature 
Casimir energy using Matsubara formalism. More precisely, the Casimir free energy can be expressed in the form of 
TGTG formula by changing the integration over £ in the formula for the zero temperature Casimir energy (|15|) to 
summation over Matsubara frequencies £ p ~ i^npksT /h. More precisely, the Casimir free energy is given by 

00 

Seas =fc B T^'Tr In (I - M(tf p )) , (23) 

J9=0 

where the prime ' over the summation means that the term with p — comes with a weight 1/2. 

The main ingredients of the TGTG formulas are the T, T matrices and the translation matrices U, V, W. Notice 
that the T or T matrix of an object does not depend on the other objects. It is derived solely from the boundary 
conditions imposed on that single object. If we consider the Casimir interaction between the object and an object 
outside it, we need to find the T-matrix of the object. If we consider the Casimir interaction between the object and 
an object inside it, we need to find the T-matrix of the object. In the following sections, we will find the T and T 
matrices of planes, spheres and cylinders under various boundary conditions. We will then discuss the translation 
matrices U, V, W between different pairs of objects. 

In this section, we focus on Casimir interaction between two objects. As in (33ll36j. one can in fact consider Casimir 
interaction between multiple objects. In Appendix \K\ we show how the approach in this section can be generalized 
to more than two objects. In that case, the Casimir interaction can still be expressed in the form (|15p but the matrix 
M has a different representation. However, the building blocks of the matrix M are still the T,T matrices and the 
translation matrices. 



III. THE T MATRIX OF A PLANE 



Choosing e z as the distinct direction, plane wave basis are parametrized by = (k x , k y ) G M 2 , with 



9?£ ut (x k) — e ik ^ x + ik yVT i \/k 2 -k 2 L z 
TE rcg 1 rcg 

A. ^ out (x,fc) =— V x ^™ t (x,fc)e z 
k 1 



k 1 



TM rog 1 reg 

A w± ' out (x, k) =^V x V x pg*(x, k)e z 

=-^-e lk - x+lk yy^ k2 - k ^ z (±k x ^Jk 2 - k\e x ± k v yjk 2 - k 2 ± e y + k 2 ± , 

Here k± = ^Jk 2 + k 2 . . 

Assume that a plane of dimension H x H is located at z = Q. In the following, when we consider the Casimir effect 
between a sphere and a plane or a cylinder and a plane, we always regard the sphere and the cylinder as inside the 
plane. Hence, we only compute the T-matrix for a plane. Let a scalar field be represented as 
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for z < (inside), and by 



poo poo ji pOO 77 

//•' / duj 2Z ^B^<(x,fc)e-^, (25) 



for z > (outside). First consider the boundary condition uip + vd z <p = on the plane. In this case, the field outside 
does not contribute to T. The boundary condition gives 



(u - iyjk 2 - k\v \ + b k± (u + iy/k 2 -k? x v\ = 0. 



where fc is given by (|S|). Passing to the imaginary frequency uj = i£, we find that 

f^m = -— = — ~ v ^ k2 + k ^ 

^' U + Vy/K 2 + kl' 

where 



£ 2 m 2 C 2 



« = V^ + — ■ (26) 



For a semitransparent plane, the boundary conditions ([3]) give 



i^/fc 2 - fc 2 + iJ fc 2 - fc 2 (a ki - 6 ki ) = AB kj 
Eliminating _B ki and passing to imaginary frequency, we find that 

T k ^£) 



b k ^ A + 2 v / K 2 + fc 2 _' 



with k given by (|26p . 

For an electromagnetic field, let 



A = H 2 l°° du> p ^ r ^ (a k - A™< ^( X , fc) + fo k -A™< °*( Xj fc) + c ki A™' ^ fc) + <i ki A™' ° ut (x, fc)) e 

(27) 



for z < (inside), and 

a=h 2 f°° duj r r My ( B ^ a te, ou t(X; fc) + ^^TM, ou t(Xj e _ iwi (28) 

»/ — oo J — oo 27T J_ oc 27T V J 

for z > (outside). Assume that the inside of the plane is filled with material with permittivity and permeability /it, 
and the outside is filled with material with permittivity e e and permeability [i e . Then the continuities of eE n . E|| , B„ 
and /z _1 B|i give 



a 



ki '- b k± = S k 



^± /„k, .k, \ _ 1 



_6 k -)=--5 k -Jfc 2 -fci, 



Vtf k j f^x _ jk^ = pk^ V k e k l 



where 



fce =V£eHeU = fl e — , 



fc 4 =j£UM.u — — 
c 



c (29) 
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Eliminating _B ki and D k± gives 



Hence 



is a diagonal matrix with elements 



where 





H\J k l 


-kl-^y/kf-kl 


5kT ~ " 




-k'i+lle-s/k'i-k'i 


c k_L 




- fc^ £ e yjk'^ k"j_ 




e iy/ k e ' 


- 4" £e\/k'i — k^j_ 




nrkx I 

t -y 


f™ \ 










+ +H e y/Ki + k'i 






+ k 2 ± -ee^n'f + k'i 




+ k'{ +s e y / Kj + k]_ 




= n e {i£)- 
c 


, Ki = rii(i^)-. 

c 



(30) 



(31) 



(32) 



IV. THE T, T MATRICES OF A SPHERE 

In spherical coordinates x = r sin 9 cos <f>,y = r sin 9 sin <fi, z = r cos the regular waves and outgoing waves solutions 
of the equation ^ are 

^ g (x, fc) = C?*ji(kr)Y lm (0, <f>), ^f(x, fc) = ^(krW^ 0). (33) 
where Z = 0,1,2,..., —l<m<l, 



are the spherical Bessel functions, and Yi m {6, <j>) are spherical harmonics given by 



Here Pp l (z) are the associated Legendre functions. The constants C[ eg and C° ut are chosen so that 

For the vector wave equation ([7|), choosing e r as the distinct direction, the TE and TM waves are 
A™' *(x,fc) =^=y V >< ^L(x,fc)re r 



C l t*(i \ ( im v (a a\ dY lm{», - i 



v /j(TTT) 1 V ' \sin9 " ' T ' 39 

IJL IV /.' i = — 

dY lm (9,<t>) im 



A i m '*( x , fc ) = , ; x V x pf ro (x,fc)re 



<? (^E2/r(*Dn.(«, *k + 7 =H into*)) 



e + -r-Jfi 



d9 sin 



Im \ 
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where I = 1, 2, . . . , — I < m < I. Here * = reg or out, with f[° s (z) = ji(z) and f° ut {z) = h^\z). 

Now we consider the T, T matrices for various scenarios. Consider a sphere of radius R centered at the origin. To 
compute the T-matrix, let a scalar field be represented as 



/oo oo j 



inside the sphere, and as 

oo I 



if : 



^ E E (^taffr fc ) + & ' m ^T(x, *)) (34) 

-OO j r> J 



Z=0 m=-Z 

outside the sphere. To compute the T-matrix, a scalar field is represented as (|34|) inside the sphere, and as 

oo / 



f = 



/oo ^ <■ 
dw E E ^™^(x,fc)e 
-OO ; ^ , 



(=0 m=-l 



outside the sphere. 

If Dirichlet, Neumann or Robin boundary condition is imposed on the sphere, and we are considering the Casimir 
interaction between this sphere and other objects outside the sphere, the field inside the sphere does not contribute 
to the Casimir interaction energy. With the boundary condition mp + vrd r ip = on the sphere, we find that 

a lm Cl cs (uji(kR) + vkR][{kR)) + b lm Cj> ut (uhf ] \kR) + vkRh{ iy (kR)\ = 0. 

Passing to imaginary frequency, we have 

lira (u - I,,i(kR) + VkRI' i (kR) 

T lm (i£) = — - = 2 LL2 

a lm (u - |) K 1+ i(kR) + vkRK' 1+1 (kR)' 

where k is given by (I26[) . 

If we are considering the Casimir interaction of the sphere with other objects inside the sphere, the field outside 
the sphere does not contribute. Using the boundary condition wp + vrd r tp = 0, we find that 

(u-V) K {kR)+vkRK > {kR) 

T lm (iO=(T lm ) i0 = 7 iyj ~, =r~ J \ m ■ 

(u — 1 1+ i(kR) + vkR1 1+ i {kR) 

If the sphere is semitransparent, i.e., it satisfies the boundary conditions ([3]), then for T lm , the boundary conditions 
give 

a lm C\ c& 3 i(kR) + 6' m C ; out /i; (1) (fci?) = A lm C T , e& ji{kR), 

a^C^kj^kR) + b^C^kh^'ikR) - A lm Cl e& kj'i{kR) = \A lm C\ cg ji{kR). 

Eliminating A lm , we find that 

, lm _ b i ™ _cr s xji(kRf 



rpirn 



a lm C° ut kj',(kR)hl 1] {kR) - kh { p' '(kR)ji(kR) + X]i{kR)hf ] {kR) ' 
Using the identity 



f l ( z )h\ 1 \z)~h ( l 1) '(z)j l (z) = --, 



z 

and passing to imaginary frequency, we find that 

XRL.iUR) 2 



\ + \RI l+ ,UR)K l+ AnRy 
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where k is given by (l26i) . 

For T lm , the boundary conditions give 



a lm C" s ji(kR) + b lm C^ ut h\ L> {kR) = S im C, out /i} 1J '(kR) 

B^Cf^kh^'ikR) - a lm Cl c& kj[{kR) - b lm Cf ut kh { ^' {kR) = \B lm Cf ut h^ (kR). 
As before, we find that 



XRK,,i(kR) 2 

T lm {i£) 



l + XRI l+i ( K R)K l+i ( K R) 
Now consider the electromagnetic case. For the T-matrix, an electromagnetic field is represented as 



/oo °° j 

rf-EE i Alm ^ ros (*> k ) + clmA ™' ros (x, 1 
'OO I 1 _ 7 



e 



l—l m——l 



inside the sphere, and as 

/oo °° j 

^ E E (« imA L E ' rCS (x, *) + 6' m A™' ° ut (x, fc) + c^A™' reg (x, fc) + d"»A™' out (x, fc) e-*- (35) 

outside the sphere. For the T-matrix, an electromagnetic field is represented as (|35|) inside the sphere and as 

/oo °° j 

^EE (s tol A™' out (x,fc)+ J D i " l A™' out (x,fc))e— ( 
-°° i=x m =_i 

outside the sphere. 

For a magnetodielectric sphere whose inside is filled with material with permittivity e$ and permeability /Xj, and 
whose outside is a medium with permittivity e e and permeability /i e , the continuities of eE„, E||, B„ and /i -1 B|| give 

a^C^MkeR) + b^Cr'h^ikeR) = A im C; cg jKfc^), 

lm C «gJ_ ^(fc^) + k e Rj[{k e R)) + b lm Cf nt — (hf\k e R) + k e Rh\ 1)} \k e R)) = A lm Cj cg 

f_l e \ / fjj e V / 

c ln lC rc g ^ ji{keR)+d lrn c out^ h (l) {keR) = C^C^^faR), 
f£ e K e Ki 

c lm Cr g ^(ji(k e R) + k e Rji(k e Rfj + d^C^y {hf\k e R) + keRh^'i^R)) = C lm Cl e ^(ji(hR) + hRj^hR)) ■ 
Here ki and k e are given by Jill). Eliminating and C' m , we find that 

where 



/T^TE n 



is a diagonal matrix with elements 

nTM 



tX e Ii + L(K e R) (il l+ l(KiR) + KiRI' l+ i(KiR) 


) - HiI l+ i(KiR) (hl l+ x(K e R) + K e RI' l+ x{^eR)j 


fJ, e Ki + i(K e R) \ 


f ^I l+h {KiR) + KiRI'^KiR)) 


- (J,iI l+ i(KiR) \ 


f ^K l+ i(K e R) + K e RK' l+1 _{n e R)^ 



e e I l+ i(K e R) (±I l+ i(KiR) - 


f KiRI' l+ i(KiRj\ 


- EiI l+ i(KiR) 


(i/ /+ i(K eJ R)4 


- K e RI'^i (/t e i£)^ 


e e K l+i ( Ke R) i 4 


- Kl RI' l+1 _( Ki R)) 


-EiI l+ i(KiR) ( 




V n e RK[ +h {n e R)^ 
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Here «j and n e are given by (1321) . 

For the T-matrix, we find in a similar way that T lm is a diagonal matrix with elements 



V. THE T, T MATRICES OF A CYLINDER 

In cylindrical coordinates x = pcos<p, y = psmcj), z = z, the regular waves and outgoing waves are parametrized by 
neZ and k z G R. They are given by 





(±K l+ i{ Ke R) + Ke RK[ +1 _{n e R)^ 


- fi e K l+ i(n e R) 


(±K l+ i( Ki R) + K t RK' l+ ,( Kt R)^j 


(j,iI l+ i(KiR) (^K l+ i(n e R) + K e RK' l+ i(n e R) 
eiK l+ i(KiR) ^K l+ i( Ke R) + K e RK' l+ i{K e R)) 


) ~ VeK l+ ±(K e R) 

- e e K l+ i(n e R) ( 


(±I l+ i( Ki R) + KiRI' lH (KiRfj 
±K l+ i( Kl R) + Ki RK' l+ i(KiRf) 


e l I l+ i{n i R) (^K l+ i(K e R) + K e RK' l+ ,{ Ke R)^ 


) -e e K l+ i( Ke R) | 


^I l+ i(KiR) + KiRI' l+ ,(KiR)) 



/„t(x,fc) ^J n {k^p)e m ^\ 

<pX (x, k) =c^H^(k^ P y n ^ z , 



and 



A ™'* (x ' fc)= fc^ VxVX ^ e2 



k± 
1 



^2 » : / 



fn(k± P )e p ^-/*(fe ± p)e^ + ^f*(k ±P )e^j e in <f+ ik **. 



Here 



k_\_ — \J k^ k^ , 



f?*(z) =J n (z), r n ut = H^{z). 
The constants C^ eg and C° ut are chosen so that 

C g ^(*C) = UO, c%*hW (*o = tf„(0- 

Consider a cylinder p = R of radius and length 77. For a scalar field, to find the T-matrix, let 

-00 J-00 Z7r „=_«, 

inside the cylinder, and 

<W |? E («^^I(x,fc) + ^^(x,fc)) e — * (36) 

-00 */ oo n= — oo 

outside the cylinder. For the T-matrix, let the scalar field be represented by (1361) inside the cylinder, and by 

/°° r°° rib °° 
-oo J — oo ^ „ 



n— — oo 
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outside the cylinder. 

First let us consider a scalar field with boundary condition utp + pd p ip — on the cylinder. This gives 
a"^C g (uJ n (k x R) + vk±RJ^(k ± R)) + b nk *C° ut (uH n {k x R) + vk x RH^' {k x R)) = 0. 

Therefore, 

b nk * uI n { 1 R) + v 1 RI' n { 1 R) 



T nk *(i£) 



uK n { 1 R) + v 1 RH'n , \ 1 R) 

a nk * _ uKjgR) + vjRH^jrfR) 
~ uI n { 1 R) + v 1 RI' n {iR) ' 



where 



7 = \ K + kl 



c 2 z h 2 ' 

Next consider the case of a semitransparent cylinder. For the T-matrix, the boundary conditions ([3]) imply that 
a nk *C^J n {k±R) + b nk *C° ut H(p(k ± R) = A nk *C^J n (k ± R), 

a «fc* C re gfcx j^ kx R) + b nk *C° ut k ± HW'(k ± R) - A nk *C^k^J' n {k±R) = \A nk *C r n cs J n (k ± R). 
Eliminating A nkz and passing to imaginary frequency, we find that 

b nk * \RI n { n R) 2 



a nk * 1 + \RI n (-yR)K n (~/R)' 
Similarly, we find that 

~ nkz _ \RK n { n R) 2 



l + \RI n { 1 R)K n ( 1 R)- 

The T and T matrices for electromagnetic fields are more complicated. Assume that a magnetodielectric cylinder 
has permittivity £i and permeability /i^, and it is surrounded by medium of permittivity e e and permeability /i e . For 
the T-matrix, the continuities of eE„, E|| , B„, M^By across the surface give 

c nk *C™^J n {k±, e R) + d^Cr^H^ik^R) = C nk *C^^ij n (k ±ti R), (37) 

rv e K e rZi 

a nk *C™^(k ± , e R) + b nk ^H^'(k ± , e R) + c nk *C^—^-J n (k ± , e R) + d nk *C^—^-H£\k ± , e R) 

tZpKi I p ±\j ftp K I o it 

nk (38) 
= A nk *C™zj' n (k± ti R) + C^C^—^-JJk^R), 

a nk 'C^^-J n (k x ,eR) + b nk *C™^H£Xk x , e R) = A nk *C* n e ^J n (k ± ,iR), (39) 

jl e jl e fli 

\i e k^. e R (j, e k ± , e R [i e fi e 

nk k ^ ' 

= A nk ^—^-J n (k ± ^R) + C nk *C™^J' n (kx,iR), 
Hik±,iR in 



where ki and k e are defined in (|29j) . The first two conditions come from the continuities of E||. They imply the 
continuity of B„. The last two conditions come from the continuities of [i~ Bn. They imply the continuity of eE„. 
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By eliminating A nkz and C z , we can solve b nkz and d nkz in terms of a nkz and c nkz , which can be presented in the 
form: 

hnk„\ / T TE,TE T TE,TM\ / h nk z \ 

(jnkzj — ~ 1 ^ c nfe z y — — I /jnTM.TE jiTM,TM I \^nk B J > 



where 



« TE «) =^ ( A™;^™' 2 + - n-fl n ( le R)K n i le R)I n{ll RY 



nkz y s; A nfez V nk * nkz j 2 ^ 2 R 2 c 4 



t TM,TM,.m 1 / A TE,2 A TM,1 . n 2 k 2 ^ 2 , ., ., 



*nk z 



ihlR 2 ^ 



mTE,TM/.« mTM,TE/.rt 1 ink z K e 2 2 2 
T nkJ (<) = T nk z (*0 = ~ T .,2 P 2„2 ("e ~ n i ) J «(7i#) . 

A™' 1 =HHeIn{liR)l' n (leR) ~ J e fJ.iIn(-yeR)l' n (jiR)> 



nk z 
TE. 
nfcz 
TM 

TM 



A ™' 2 =H»eIn{liR)K{ le R) - le ^K n { le R)l' n ( 1% R), 
,TM,1 



A „fc!' =HSeIn^iRK^eR)-7eeiIn(7eRK^iR), 

f: 2 =^e e I n ( lt R)K' n ( le R) - le eiK n {j e R)I^iR), 

n 2 h 2 P 2 

a TE,2 a TM,2 . « K z? / 2 2\2 7^ / r>\2 r 1 D \2 
Iile rl L 



le = \/ K l + k 2 , Ji = \Jk 2 + 



k 2 . 



where Kj and n e are given by (J32J). 

For the T-matrix, we find in a similar way that 



^TM,TM/.m 1 ( TTE,2TTM,1 . n 2 k 2 £ 2 2 2 2 2 



1 ZTlk-rKij 



A™' 1 =7eft^(7e«)<(7^) - 7i^ n (7ifl)4(Tefl), 

A™/ =^iK n (j e R)lU7iR)-7iPeIn(liR)K^eR), 

;tm,i 



A nfc.' = 7 e£i^n( 7eJ RX( 7 i^) " 1iE e K n (7iR)K( le R), 

A™' 2 = 7 e £ i {leR)l' n {liR) ~ ^U^K^R) , 

„2i,2t2 

TTE,2TTM,2 . / 2 2\2 f /• D \2 r / D \2 

= A nfe J A nfc ; + 7 2 7 2 E 2 C 4 ( n e~"i) ^«(7ei?) / n (7i-K) 



*nk- 



Notice that the T and T matrices are diagonal if and only if n$ = n e , i.e., the media inside and outside the cylinder 
have the same refractive index. In this case, 

T TE,TE = liH.eIn{liR)I' n {"1eR) - J e ^jIn(leR)In(%R) 

nk > jmJ n ( 7i R)K' n ( 7e R) - ^K n ( 7e R)I^R) ' 

j,TM,TM = li£eIn(liR)I' n (leR) - le£jIn{leR)I'n{liR) 

nk * jie e I n ( 7i R)K' n { le R) - le e t K n ( le R)I^R) ' 

^TE,TE ^ 1eVtK n (leR)K' ll (j l R) - Jjfi e K n ( 7 j -R) K' n ij e R) 
^TM,TM _7e£ 

iK n (j e R)K' n { 7i R) - ^Knij^K'^R) 
nk * 7e£iJ ft-„( 7eJ R)74( 7iJ R) - lt eJ n ( lt R)K^ le R) ' 
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VI. CASIMIR EFFECT BETWEEN TWO SPHERES 



In this section, we derive the TGTG formula for the Casimir energy between two spheres using the prescription 
in Section [TX| Since we have discussed the T and T matrices in Section IIV1 now we only have to find the translation 
matrices U, V and W. These matrices have been derived in a numbers of works [48l 453| . Since the operator approach 
used in [48[ can be used for other scenarios, we will discuss in detail the approach of |48|. 

Consider two spheres, one is centered at O = (0,0,0) with radius R\ and one is centered at O 1 — (0,0, L) with 
radius i?2- Assume that R2 > R\. 

Define a differential operator Vi m by 



2l + l (l-m)\ f d x +id y \ m M (d z 



An (I + m)\ \ ik J \ik 

Here m > 0, Pi{z) is the Legendre polynomial of degree I and P[ m \z) is its m-times derivative. It follows from the 
definition of spherical harmonics [54J that 

^ m e lk ' r = y ;m (^.,(/) fc )e lk - r . 

Here k = k x e x + k y e y + k z e z , r = xe x + ye y + ze Zl k x — k sm6 k cos 4> k ,k y = k sin^fc sin <j) k ,k z — k cos 9 k . On the 
other hand, it can be proved by induction that 

Vi m jo(kr) = i l ji{kr)Yi m (B,<l>), Vi m hg\kr) = i l (kr)Y lm (9 , 0) . (41) 

Using the integral representations 

Mkr) = S -^j^- = ^ J J d& J" d9 k sin0 fe e lkr , 



exp(ifcr) 1 r „ r „ e **.*+**»v±y* a -fcS-fcS* 



ikr 27r 



x k^/k 2 — k 2 — k 2 



where A /fc 2 — k 2 — ky = i^Jk 2 + k 2 — k 2 if k 2 + k 2 > k 2 , we have 

TTOg n27T 



£> re 5 /.^7T />7T 

^L s ( x - *>) =4^r^m y o y o sin # fee i 



1 

'47Ti ( 



2tt 

d(j>k 
Jo 



(42) 



d9 k sm9 k Y lm (9 k ,<f> k )e ik - r , 



and similarly [551 ]: 



^>out poo poo ik x x+ik y y±iy/k 2 —k2 -k?z 

^(x,fc)=^ r / dfc x / dfc,lW# fc ,0 fc ) , , z%0. (43) 

'-00 J-00 k^/k 2 — k 2 — k 2 



Here kcasdk — k z — ±y ^ 2 — k 2 — k 2 for z ^ 0. The integral representations (|42|) and (l43l) will be very useful later. 
In fact (|42[) is equivalent to the well-known formula 

00 1 

e lk - r = 4vr^ J2 ^^kr)Y lm {eA)YC m {9 k Ak)- (44) 



Now, consider the expansions which define the components of the translation matrices: 

oo I 

^,(x-L,fc)=£ £ VWw(-L«(x,fc), 

1=0 m =-l 
oo I 

tf#,(x-L,fc)=E E £W™'(-I«(x,fc), 

oo 

r=o m'=-V 

oo Z 

^(x' +L, A) = E E W Vm , M (±)<p1%,{x!, k). 

l'=0 m'=—l' 

Since 

/->reg „27T ^.7r f — IV™ 

(fl»m»^ s ) (0) y d0 fc y smff fc y Im (fffc,^ k )yi// m «(ff fc ,^) = ^-^.i"^",-^ 168 , 

apply the operator Vi- m to both sides of (|45|) and set x = 0, we find that 

VWm'(-L) = ^5f ("I)™ (n-m^/) (-L) 

= (-l) m / d0 fc / ^ fc sin0 fc y i! _ m (^,0 fc )r rm ,(^, ( /» fc )e lk - L . 
Jo Jo 

Here we have used the fact that (Pj.-mV/w) ( — M = (— (Vi t - m (p T r?!L,) (L). Now using the identities [Hf 
Yi- m {9 k , (f>k)Yi< m '(9k,4>k) 

oo Z 



V V r IV"", / (2Z + l)(2Z' + l)(2i" + l) / i I' Z" \ (I V l»\ , . 

<L V 4^ Urn m' -ml') U J y '"™"(^> 

' ' '" \ . f_iy+i'+i" i ' '" 



-m m — m 

and the fact that 

i r i" \ (i v i" 

m -ml m") ^0 

is nonzero only if I + V + I" is even, \l — l'\ < I" < I + I', and m — ml + m" — 0, we find that 

Am 1 " /(2Z + l)(2/' + l)(2/" + l) 



v lm>Vm ,(-L)=(-ir E E 



2+2' J" , 



'=|J— i'l m"=— I" 



X (m -m' to") (o ) ( L ' fc )- 



In the case L = Le z , we have 



r >;; i ,,(LJ, ) =t7;V^^^(^)^,o. 



Passing to the imaginary frequency, we find that 



2+2' 



' K l"=l«-J'l 
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where 



Hj;. m = y/(2l + l)(2l' + l)(2l" + 1) 
In the same way, (|49l) and (l4l)l) imply that 

U lm , Vm .{--L) =^(-l) m (Vi,- m ^) (-L) 



I V l"\ I V V 



m -m 



(51) 



7r(-l) ( ' +1 (-l) m / dk x I dkyY^tfk^Yvrtlptofa)- 



oo poo 



— oo J — oo 



ik 



.x+ik y y+i^k 2 — k1 — k^z 



Jc Jc 2 — h 2 — Jc 2 
nj \ I t\j x y 



l+V l" 

=(-i) (,+i (-i) m e e (-ir'X 

j"=|2-J'| m"=-l" 



no 



2tt¥ /(2Z + 1)(2Z' + 1)(2/" + 1) 



4tt 



Z Z' Z" \ I V I 



to —to to 



" (o Wrm"(^k). 



Specializing to L = Le z and passing to imaginary frequency, we have 



U lm j, m ,(~Le z ,tO = (-l) l ' +m 6 m , m , x /^ z E H i ;;, vm K lll+h {nL) 



i+i' 



Analogously, (gSJ) and flU) imply that 



i"=\i-v\ 



2+2' 



2" = |2-2' 



Finally, notice that for general L, 

oo V 



('=0 m' = -2' 

=E E (-!)" +m 'E E (- ir w ^/ (2I+1)(2i/+1)(Mff+1) 



i'=0 m' = -V 



47T 



2»=0 m" = -f" 

.to -to' -to") (o ) ^"" l '' (L ' fc) ^'™' (x '' fc) - 



X 



Interchanging x' and L, we have 



oo I 



oo 2 



■•Pi, 



t (x'+L)=E E (-ir'E E (-i) 1 " 4 *"^ 



2'=0 m'=-2' 



2"=0 m"=-l" 



(2Z + 1)(2Z' + 1)(2Z"+1) /Z Z' Z" \ /j Z' Z" 



4tt 



to -to' -to" / 10 



x^ e 4„(L,fc)^(x',fc). 
Compare to (|48|) and specialize to L = Le z , we find that 



I i+i' 

W Vm , M (Le z ,it) = (-l) m S m , m J JjL (-l) l "H l u l{m I lH+k (KL). 

" 2" = |2-2' 

Hence, for a scalar interaction between two spheres, the Casimir interaction is given by 



E Cas = — d£Trln(l-M(i£)), 
27T Jo 
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where the trace Tr is 



£> 

m=— oo l=\m\ 



m=—c 

and M is a matrix diagonal in m. When the two spheres are outside each other, 



(52) 



M,, 



OC 

/— |m| 



7T 

2kZ 



i"=|i-f| 



7T 

2kI 



When sphere 1 is inside sphere 2, 



E 

i=|m| 



2kL ^ 

i" = |Z-i 



TV 



I'+l 

E < m ^ + i(^)- 

f"=U'-f| 



. V ff'',' Ij„,i(kL). 

f" = |Z'-f| 



Here the fact that H l ~ and i/f r are nonzero only if I + I + I" and I' + I + I" are even have been used to write 

t/;m u l\m 

(— 1)' +l as (— 1) + , and the cyclic property of the trace has been used to get rid of this sign factor, k is related to 
£ by (|26p . T\ m is the /m-diagonal component of the T-matrix for sphere 1. T^" 1 and T% are components of the T 
and T matrices for sphere 2. They have been derived under various boundary conditions in Section HVl 

The exact Casimir interaction energy for two spheres that are outside each other has been considered in a number 
of works. For two Dirichlet or Neumann spheres, it has been considered in [23, IH, H3- For two Robin spheres, it 
has been considered in [33l. [34l |57|. For two semitransparent spheres, it has been considered in [46|, S3]- The Casimir 
interaction energy for two spheres with Robin boundary conditions, where one is inside the other, has been considered 
in [57j • Besides these scenarios, there are other interesting scenarios whose exact Casimir interaction energy can be 
obtained from above: 

• Two semitransparent spheres, where one is inside the other. 

• A Dirichlet/Neumann/Robin sphere and a semitransparent sphere, where the two spheres are outside each other 
or one is inside the other. 



Next, we consider electromagnetic fields. Let C be the operator 

1 



-r x V 



smt 



and define 



so that 



--CY lr 



sm( 



dY, 



Of) 



A™'*(x,fc)=Cf/r(*r)X, m ( 



Let 




{CPlra ~ VlmC) 



- m)(l + m+ rjVi, m +l + y/(l + m)(l-m+l)Ti, m -i 

— ime z Vi 



(53) 



It follows that 



-m)(l + m + l)Vi, m+ i - y / {l + m)(l-m+l)VL m -i 



A™' rcg (x, k) = Cr^ l T lmJO (kr), A™' out = CrH- l V lm h {kr) 



(54) 



In [48j], it was shown that 



This implies that 



A™' rog (x,fc) =^— l j d<f> k J d8 k sm8 k X lm (6 k ,0 k )e lk - r , 

-, ^reg „2tt rix i 

A™' rcg (x, k) =±V x A™' rog (x, k) = g-^ J d<f> k J d9 k sin# fc - x X lm {9 k , <t> k )e< 



A™' ° Ut ( x ) fc ) t / *x / dk y Xi m (9 k ,(f>k)- , = , z^O 



/-"out poo roc ikxX+ik y y±iy/k 2 — k 2 —k 2 z 



2iri l 



>- k^/k 2 — k 2 — ky 



Qout roc roc ^ ^ik m x+ik y y±i^J k 2 ~k 2 —k 2 z 

A™' ° Ut ( x i fc ) =7T — 7 / *x / dk y — x Xi m (9 k ,(j) k ) = , 

" - '-oo fc k^/k 2 — k 2 — k 2 



Now consider the expansions that define the translation matrices 

l 



A™' -(x - L, *) = £ E (O™ (-L)AL E ' -(x, lb) + 0^(-L)AL M » ^(x, *) 

/— 1 rn— — l 
oo I 

A™; " cs (x - L, fc) = 2 E (0™(- L ) A ™- ^ fc ) + C% M (-L)A™< rcs (x, fc) 



Z=l m=—l 

Using the relation ©, it is easy to see that 



t^TE.TM _ ttTM.TE y-TM.TM _ ^TE,TE 
Irtijl'm' hn,l'm' 1 Ira^l'm' lm,l'm'' 



Since 



■ v , m „ ■ A™' rog J (0) =^— l j d<j> k / dO k sm0 fe X r , ro „(0 fe ,<^) ■ X lm {9 k ,4>k) = ^T^M'"^ 
{V v , m „ ■ A™< rcg J (0) =-±-r l J dcp k j d9 k S in0 fe X w (0 fc ,cfe) • ( — x X iro (0 fe ,<fe)J = 0, 



applying the operator "Pz.- m - to both sides of (1581) and setting x = 0, we find that 

0™(- L ) ("D m (n- m • A™-) (-L) 

=(-l) m / d0 fc / d^sin^X^^fe^^-X^^,^)^^, 
Jo Jo 

(-l) m (TV*. • A™; reg ) (-L) 
=(-l) m+1 #ft jT <W fc sin0 fe X;,_ ro (0 fe) fc ) • ^ x X,, m ,(0*, <f> k ) ) ( ' k i 
Here we have used the fact that 

(TV™ • A™-) (-L) = • A?5 rCS ) (L) 

and 

(Pi,-m • A™; rcg ) (-L) = ( n _ ro . A ™;-) (L). 
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In Appendix [B] we verify that 



x^ m{ e M • WM*) - E ^ 2y /i(i + m i> + i) M) 



and if L = Le z , 

r 2ir 



acp k 
o 



,/ (2I + l)(2i' + l)(2i" + l) f| V 
x V 4^ U ~m' m" U J m " (0fe ' ^ fe) ' 



d0 fc sin0 fc X,,_ m (0 fc ,0 fc ) ■ I ^ x Xi w (0 fc ,fo) ) e lkL 



7 T / > 27T />7T 

vmkL 



k L 



^(/ + i)z'(z' + i)y V 

Compare to the scalar case, we find that 

l+V 



fc / d8 k sin6 k Yi, m ,(9 k ,(j) k )Yi- m (9 k ,(j) k )e l 



(«L), 



where 



In the same way, we find that 

/+/' 



/// i/ // 1 ^ H ]tfi n' / 



(61) 



(62) 



A£ = ^ + + + A«/. m = (63) 

2^/Z(Z + !)/'(/' + 1) ' v^+IFcF+i) 



w, mVm (Le z ,io=(-ir +1 wy^ E t- 1 )'" (a^L A A{f)^ J '"+i^- 

Gathering the results, we find that for electromagnetic interaction between two spheres, the Casimir interaction 
energy is given by 

Ecm = 2^J ^ Trln ( 1 - M ( i 0) 5 

where the trace Tr is 

oo oo 

Tr = E E tr > ( 64 ) 

m=— oo Z=max{l,|m|} 

and the trace tr is a trace over 2x2 matrices. M is a matrix diagonal in m. When the two spheres are outside each 
other. 



i+i 



A'r -A, 



r ri i n V r,/ i ; ?! V H;m |J / 

!=max{l,|m|} = 1\ \ ' " / 



I'+Z 



2 V 2kL ^ \A„ T A 1 - / " ;m +5 

r«_ i»/_ri \ vi\m i'i / 



i"=|i'-i| 
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When sphere 1 is inside sphere 2, 



l+l 



r=max{l,|m|} i"=|Z— 1\ \ " ,m " ' 



2 (5 ****** 



2kL ^ \A,,r A' - 

Here k = n m (i£)£/c, where n m is the refractive index of the medium between the two spheres. The matrices T 1 ™ for 
sphere 1 and T l ™ , T l ™ for sphere 2 have been derived in Section IIVI 

As in the scalar case, the electromagnetic Casimir interaction between two spheres has been considered in a number 
of works. For two perfectly conducting spheres that are outside each other, it has been considered in [32l. l34l |45|. l57l. [58jj . 
For two magnetodielectric spheres that are outside each other, it has been considered in_|32|, [3J, [58( . For two perfectly 
conducting spheres, where one is inside the other, it has been considered in [57], [H, |60|. For two magnetodielectric 
spheres, where one is inside the other, it has been considered in [60| . 

VII. CASIMIR EFFECT BETWEEN A SPHERE AND A PLANE 

In this section, we consider the Casimir effect between a sphere and a plane. This scenario has been extensively 
studied for its experimental value. It is more appropriately regarded as one object (the sphere) is inside the other 
(the plane). Assume that the sphere is given by x 2 + y 2 + z 2 = R 2 , with center at O — (0,0,0) and radius R; the 
plane is given by —H/2 < x, y < H/2, z = L, with center at O' = (0, 0, L). 

Since the T, T matrices for a sphere and for a plane have been considered respectively in Section IIVI and Section 
IIIII let us consider now the translation matrices V and W. In the scalar case, they are defined by 

^°f(x-L,fc)=£ WJ-L)¥C g (x,fc), (65) 

1=0 m=-l 

^(x' + L,fc) = H 2 / p ^W kliim (L)^(x',fc). (66) 
J-oo 2vr J_ oa 2tt 

For Vi m ,k ± (— L), changing z to — z in (J3U), and multiplying by e ' l \/ k2 ^ k2 L L , we find that 

oo I 
1=0 m=-l 

Compare to (|65|) . we have 

V lmM± (-Le a ) =47r(-l) m ^„(^,</) fc )e J V^ii 

For Wk Xl ; m (L), compare (|66l) to (H31 . we obtain immediately that 

^ l 2l + l(l-m)\ ( yfW^WA ( k x + ik y \ m iy ^zkT L 

H 2 ky/k 2 - k\ V ^ (l + m)\ 1 \ k )\ k ± J 

Consequently, we find that the scalar Casimir interaction energy between a sphere and a plane is 

E c ,s = tt / d£Trln(l-M(zO), 

27T Jn 
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where the trace is given by (|52p , and 



oo poo 



— oo 



M lrn , Vm ,(u) = - ttT 1 " 1 I dk x I dk . yi,-m(0fc,fo)ri'm'(0fc, <j) k )T^e 2 ^ k2 - k ^ L 



k^/k 2 - k'j 

^+1 Tjl pOO pOO i 

= -ir5 mm ,T lm Y - " V " / dk x / rffc v - ^» o (0 fc , ( /) fc )T k ^e 2 V fc2 -^ L . 



l"=\l-l 



Here ([50)1 and (|5T|) have been applied. The T-matrix component T lm for the sphere has been derived in Section ITVl 
and the T-matrix component T k± for the plane has been derived in Section Hill If the plane is subjected to Dirichlct 
or Neumann boundary conditions, T ki = 1 or — 1, in this case, we obtain from (|43p that 



2+2' 



M lmJ , m ,(iO = ±5 m , m 'T lm J-?- ]T H l u,. m K lll+ i(2 K L). 



4kL 

i"=\l-l>\ 

In the general case where the plane is semitransparent or is subjected to Robin boundary conditions, 

l+V "°° '- / ^K 2 + k 2 ± 



2 v ,^_ n Jo n^FTkJ \ J 



Making a change of variables 

k±_ = k sinh 

we find that 

2+2 -.oo 

l"=\l-l'\ Jo 
= ^6 nhm ,T lm J2 dzP v ,{z)f^e- 2KLz . 

For a Robin plane, 

vl — k cosh u — KZ 



2kL cosh 



u + k cosh 9 u + kz' 
whereas for a semitransparent plane, 

f ^ (if) = - = - 

v ; A + 2kcosIi# X + 2kz' 

In the formulas above, K is related to £ by (|26|) . 

The exact Casimir interaction between a Dirichlet sphere and a Dirichlet plane has been considered in [27], HH, H3, 
l6ll EH . The interaction between a Robin sphere and a Dirichlet plane has been considered in [(33| . Besides these 
scenarios, there are other scenarios whose exact Casimir interaction energy can be obtained from above: 

• A Dirichlet/Neumann/Robin sphere in front of a Robin plane. 

• A semitransparent sphere in front of a semitransparent plane. 

• A Dirichlct/Neumann/Robin sphere in front of a semitransparent plane. 

• A semitransparent sphere in front of a Dirichlet/Neumann/Robin plane. 

Note that the formula for the Casimir interaction energy between a Robin sphere and a Robin plane is considerably 
more complicated than the corresponding formula for a Robin sphere in front of a Dirichlct/Neumann plane. 
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For the electromagnetic case, the translation matrices V and W are defined by 

OG / 



A™' -(x - L, k) £ (^™(-L)A^ rcg (x, *) + C^ B (-L)A™' rcg (x, *)) , 

/— 1 rn — — l 

oo / 

A™' «*(x - L, fc) = 2 E (^™(" L ) A ™- ^ fc ) + ^™(-L)A™' rCS (x, *)) , 



l—l m—- 



\ TE, out / /.f- j \ rr 2 f°° dk x f°° dk y / tE.TE/t \ A TE, out , i ,\ . W TM,TE/ T \ . TM, out/ / , \ 

A im (x+L,k)=H J —J -± [W k J m (L)A k± ' {*,k) + W k J m (L)A kx - (x , fc) 

A TM. out/ / . t t \ rr2 Z" 00 ^kt i^" 00 ^^J/ / jrrTE.TM/r \ a TE, out/ i ,„TM,TM /r w TM, out/ / , \ 

A lm [X+L,k)=H y (^k^m ( L ) A k/ ( X > fc ) + W kx,Im ( L ) A k^ ( X ' fc ) 



As in Section fVTl the relation ([9]) implies that for Z = V or W, 



? TE,TM V TM,TE 17TM.TM rrTE.TE 



For V, (H7J) and O imply that 

C£ E (- L ) =|S (n- • A ™ ,res ) (-l) 

<k! M (-L) =|ei(-l)" 1 (P,- m • A™-) (-L) 



— -B x -' 



ik x x+tk yV -iy/k 2 -k 2 ±z l k % \A 2 - k ± , kyy^k 2 - fej fc x 



fcfc I 



fcfc I 



Using the fact that 



(67) 



(68) 



r=-L 



tti kki k 



-ea. 



where cos 



s/W^k 



and the relation (|55|) with fc 2 = —\Jk 2 — k^_, we find that 



Ck T V^ 2 ) =4 7 r(-l)<+™X i> _ m (0 fe ,^) • (-fe^JeV* 5 ^ 



4ttz dYi 



j - m {0k,4>k) l y /k 2 ~k 2 1 



,/TE,TM/ r \ 



=4 7 r(-l) ;+m X i ,_ m (0 fc ,</> fe ) • (- e?fc ) eV 

4?r * m ,/ u/k 2 -k 2 L 



k 2 -k 2 L 



+ l) sin(9 fe 



Here 



and we have used the fact that P, m {cos9 k ) = (-l)' +m P, m (cos6' fc ). 
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For the matrix W, we obtain from (15 70 that 



TE, out^„/ 



(x' + L,fc)=- 



1 



4 v /7(zTi) 7- 



oo />oo 

dkx I dky 



oo «/ — oo 



-Yi m {O k ,(j) k )ee k e 0fc 



sm w fe 

ik x x' ' -\-ik y y' ' -\-iyJ k 2 —k? L (z > +L) 



(69) 



Since 



e e fc 



k^ "\J k^ fc^ fc^ i 



j Z k*y ^ kj^ 



kk i 



compare to (f6"5)) . it follows immediately that 



„ 7 TE,TE rro v i 7T 2 dY im (0 k , 0fc ) j^pZfcsr ' L 



TjrTM,TE/ r X 



i? 2 V^+T) V fc2 - fc 2 sin <9 fc 
Hence, the electromagnetic Casimir interaction energy between a sphere and a plane is 

ft 



Cas 



2vr 7 



d£Trln(l -M(i£)) , 



where the trace Tr is given by (|64p . and 



I in 



y/l(l + l)l'(l' + l) 



oo poc 



1 



dkrj* I dky — 
-oo J —oo k-\J k k 



( dYi ^ m (e k ,(f> k ) 

d9 k 

-Yi^ m {e k ^ k ) 



-Yl^- rn (0 k ,4> k ) 



sm V k 



x T k 



( d Yl , m ,(e k , M JaL Yl , m , {eM 

o9 k sm 6k 
ml . . dYi> m t(6 k ,4> k ) 

\^0- k Yl ' m ' {9k ^ k) 0fc / 



sin fc 

dYi^ m (e k ,(j>k) 
d9 k 



\ 



„2iy/k 2 -k 2 ± L 



Here fc = n m uj/c, where n m is the refractive index of the medium between the sphere and the plane. The T- matrix 
component T lm for the sphere has been derived in Section IIVI and the T-matrix component T kj - for the plane has 
been derived in Section Unl If the plane is perfectly conducting, we take e e — > oo in (f5T|) . which gives 



1 

-1 



Then 



where 



Wm'M = vrT' m / dk x / dfe 



OO P oo 



oo «/ — oo 



'k^/k 2 -k'i \- A 2 -A X 



Ax A 2 \ 2i ^5T 



Ax 



A, 



1 



dYi^ m (9 k ,4> k ) dYi, m ,(6 k ,(f> k ) 



d6 k 



+ 1)1' (V + 1) Vsin 2 fe 
y/l(l + l)l'(l' + l) \sm6 k l ' fe ' w ddk ^ S in9 k lm ^ k ^ k > 



d9 k 

dYi- m (9 k ,4> k ) 
dO k 



-^~l,-m{6k,4>k) ■ ( -T x Xj/ m '(0fc,<Afc) 
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Compare to the translation matrix \Ji m ,i'm' (— L) for the case of two spheres, we deduce that when the plane is perfectly 
conducting, 



' i"=|j-j'| ^ u '' m lv ' 

where A l w is given by (joUj) . and 



y/l(l + + 1) 

Here k = k to = n m (i£)£/c. In the general case, Mi mi j/ m '(w) is nonzero only if m = ml . Passing to imaginary frequency 
and making the change of variables k± = ft sinh we have 



IT (ift-X t}¥l teRM) (l-m)\(l'-m)\ lm f°° a -2n m L cosh 9 

Wm'W)-*m,m — \/^ + W + 1) (/ + m)!(/' + m)| T X 

777 \ / / 77?' / 

sinh 0P, m ' (cosh 0) P m (cosh0)\ / sinh0P,™ '(cosh0) — — PJJ 1 (cosh0) 

sinh0 I ^pkj_ sinh0 

~^Trfl Pim(c0sh ^ sinh 0P/"' (cosh 0)/ I _ILR"'( C0S y) sinh 0P™'' (cosh 0) y 
The matrix T kj - is a diagonal matrix. If z > L is a dielectric half-space with permittivity e e and permeability /i e , 



Mm \ 


/f ■ 


f sinh 2 6 


' — /i e cosh 




/ n% 

1 ™m ' 


f sinh 2 6 


' + /i e cosh0 


£m \ 




f- sinh 2 


— e e cosh 






f sinh 2 6» 


+ e e cosh0 
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The exact electromagnetic Casimir interaction between a sphere and a plane has been considered in a number of 
works. For perfectly conducting sphere in front of a perfectly conducting plane, it has been considered in [35l. |62M65| . 
For a dielectric sphere in front of a perfectly conducting plane, it has been considered in [35[ . For a magnetodielectric 
sphere in front of a magnetodielectric half-space, it has been discussed in (3(| HI, H3] • 



VIII. CASIMIR EFFECT BETWEEN TWO PARALLEL CYLINDERS 



In this section, we derive the TGTG formula for the Casimir interaction between two parallel cylinders of length 
H . Assume that the radius of the cylinders are R\ and R2 respectively, and the centers are at O = (0, 0, 0) and 
O' = (L, 0, 0) respectively. The symmetric axes of the cylinders are parallel to the z-axis. 

The T and T matrices of a cylinder have been derived in Section |Vl Now we consider the translation matrices 
between the cylinders. These have been derived using different approaches, and they are essentially the addition 
formula for Bessel functions. Here we mimic the approach for spheres. For our application, we only consider a 
translation by a vector L that is perpendicular to the z-axis. It is obvious that all the translation matrices should 
then be diagonal in k z . 

Consider the well-known formula: 

00 

e ik*x+iky V = ^ i"J n (fc ± p)e in *e- m ^, 

n— — 00 

where 

x = pcoscf), y = psin</>, 

k x — k± cos 4>k, k y = k± sin <fik . 



From this formula, we can deduce that 



/-jreg /»2tt 

V™1 (x, k) = CZ*J n (k x p)j n *j k " = / # fc e in * fc e ik ^ p+i ***. 



Here = + fc^e,,, p = xe x + ye y and fcj^ = \Jk 2 — k 2 = yk 2 + k 2 ,. 
Define the operator Q n by 

d x +id y \ n 



Qn = 
Q-n = 

where n > 0. It is obvious that for all n, 



ik± 
d x - id y 
ik i 



Hence, 



nrcg r2n 



This identity can also be proved by induction. In fact, using induction, one can also prove that 

r n Q n H^(k ± p) = H^(k ±P )e m ^. 

On the other hand, one can show that 



I f 00 ±iyj k 2 ± -klx+ik y y 

H { 1] (k ± p) = - dk y , x^O. 

'k 2 , -kl 



7T 



V 



Hence, 



/iout poo ±i^/k^^k^x+ik y y+ik z z 

<fX (x, k) = ^ dk y e m ^ =_ , x % 0. 

Ik 2 - k 2 y 



Now consider the expansions that define the translation matrices: 

oo 

^ C ,i>-L,fc)= ]T K,„K-LKX(x,fc), 

n— — oo 
oo 

vC£. 2 (x-L,fc)= ^ ^(-IO^x,*), 

71— ■ — OO 
OO 

( /? °t(x' + L,fc)= ]T C/„,,„(L)^(x',fc), 

n' — — oo 
oo 

^£(x' + L,fc) = W n ,, n (L)<p<$l(x',k), 

n' — — oo 

with translation vector L = L x e x + L y e y . Since 

(Q„»C fe g J (°) = / tyke*"**?"** = cz*i-*8 n „,-n, 

applying the operator Q_„ to both sides of ([74]) . we find that 

7™ 

K, n ,(-L)=^(fi_„^J(-L) 

1 /-27T 

=_ / d0 fc e l (™ -«)^e ik ' L 



2tt .,. 

=(-l)"'->^(L.fc)- 
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Here we have used the fact that (Q_ n <^™,| ) (— L) = (—1)™+" {Q- n ip T ^ k ) (L). Specialize to L = Le x and passing to 
imaginary frequency, we find that 

V n , n .{-Le x ,iO = (-l) n '-"/„_ n ,( 7 i), 

where 

7 = Vk 2 + fcf, 



and fc = iK. In the same way, we deduce from (|78[). (|75|) and (|76| that 

[/^(-Le^) =(-l) n, « n _ n ,( 7 L) ) 
U n ,, n (Le x ,i£) ={-lf Kn^.^L). 

In fact, for the latter, we have 



¥>3£ (x' + L, k) = (-1)"'^,^ (L, (x, fc). 

n' — — oo 



Interchanging x' and L gives 



, „out 



;t(x' + L,fc)= (-l)"'^ U -„'^(x',fc)^(L,fc) 

77.' = — CO 
OO 

Compare to f[77|. we obtain 



n — — oo 



W^Le*^) = (-1)""" ^ n ,(Le x ,iK) = (-l) n " n W^)- 
Consequently, for a scalar interaction between two cylinders, the Casimir interaction is given by 

fc TT pOO pOO 77 

E Cas = TT / / ^Tr In (1 - M(»0) , (79) 

where the trace Tr is 

oo 

Tr= E • (80) 

71— — OO 

We have used the fact that the T,T and the translation matrices are diagonal in k z to take out the trace over k z . 
When the two cylinders are outside each other, 

oo 

M n>n ,(iO =T x nfe * J2 K^^T^K^^L). (81) 

n— — oo 

When cylinder 1 is inside cylinder 2, 

oo 

M n , n ,(iO =Tf* E In-nhL)f^I^ n ,(jL). (82) 

h— — oo 

Here the cyclic property of the trace has been used to get rid of the sign factor (— l)" - ™'; 7 = ^ n 2 + k 2 and k is 
related to £ by (|2l))) . T" is a component of the T- matrix for cylinder 1. and T£ are components of the T and T 
matrices for cylinder 2. They have been derived under various boundary conditions in Section IVl 



Next, we consider the electromagnetic case. Notice that 



A™>*(x, k) =C* n (J±ft(kj_p)e p /*'(fc±p)e^ e^+^ 

=C;^ (£ + i(fciP)e i(n+1) ** +^-i(fe±p)e <(n - 1) **) 
+ f (£ + i(*±p)e i(n+1) ** - /'-i^p)^"- 1 **)^ 
where * = reg or out, f^{z) = J n (z) and f° ut (z) = h£\z). Define 

Qn = y - Qn-l) + J (Qn+1 + Qn-l) ■ 



ik z z 



Then 



Since 



A™;*(x,fc) = r"c;s„/ *(^ P )e 



ik z z 



ikj_ p _ 



(i sin^e^ - icos ^e^) e m0fc e lki ' p , 



+ !)0fc _|_ e i(n-l)4>k 



we have 



p-\-ik z z 



2ni 

Qreg 



6^ (i sin 4>kG x — i cos 0fce y ) e m<pfc e 



TM, reg/ 



v ,. Hx,fc)=^VxA™; res (x,fc) 



/ # fe — x (tuin^e* -icosc/) k e v )e m < l,k e lk± - p+lk * z 
Jo k 



Similarly, 



TE, out 



TM, out 



(x.fc) 



(x,fc)=— 



dk y (i sin ^>fee x — i cos 0fce y ) e m< ^ fc 



a ±i-^/fcJ -k^x+ik y y+ik z z 



dk y [~ k 



k z cos 0fc k z sin 0^ fc_L 



k ~ Gy + T Gz 1 e 



Consider the expansion that defines the translation matrix V: 

oo 

A™f g (x-L,fc) = J2 C; TE (-L)A™ : -(x,fc) + <„ M ^ E (-L)A™^(x,fc), 

n— — oo 

oo 

A^; rcs (- - l. *) = E C E '™(- L ) A ™; rcg (x, fc) + V™'™(-L)A™'"*(x, fc). 
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ti r -ii r Jnl n , ,rTM,TE T ,TE,TM , T/ TM,TM T rTE,TE AT 

It follows from ([9]) that V: = V: and V n _,' = „; . Now, since 



Qr 



A nkl 1CS ) (°) =7T-^ / <#ft (isin<f>ke x - icas<f> k e y ) ■ (isin0 fe e x - icos(f> k e y )e 

I ZTTI Jr. 



i(n"+n)0fc 



Qn" • A™-) (0) =gl j( V (isin^e, - ico sM . (_^^ e , - ^^e, + %e») e^"+^ = 0, 

(85) 



we find that 



^TE,TE 



Hence, 



V n ,n'(-Le x ) =V n ,n' (~Le x )I, 
where I is the 2x2 identity matrix. In the same way, we find that 

U n ,n>(-Le x ) =U niV ,'(-Le x )I, 
V n ' tn (Le x ) =U n > tn (Le x )I, 
W n ,. n (Le x ) =W n ,, n (Le x )L 

Consequently, the electromagnetic Casimir interaction energy is similar to the scalar case given by ([29|, ([81]) and (|82|l . 
except that we have to replace T™ fe * , T 2 " fc * and T" fe * by 2 x 2 matrices Tf z ,Tf J and T™ fc * derived in Section El and 
7 is replaced by j m = yj + k 2 , with K m = ra m (i£)£/c. n m is the refractive index of the medium between the two 
cylinders. 

The exact Casimir interaction energy between two eccentric Dirichlet/Neumann/perfectly conducting cylinders have 
been derived in [4l] - |43j using the mode summation approach and in [36| using the multiple scattering approach. For 
two Dirichlet/Neumann/perfectly conducting cylinders that are outside each other, the Casimir interaction energy 
has been considered in p], |3l], S, [H, For two magnetodielectric cylinders, the Casimir interaction has been 
considered in [70|. For two semitransparent cylinders that are outside each other, the Casimir interaction energy has 
been derived in |46l. l47j . In addition to these scenarios, the results in this section also provide the exact formula for 
the Casimir interaction energy in the following scenarios: 

• Two Robin cylinders that are outside each other, or one is inside the other. 

• One Dirichlet/Neumann/Robin cylinder and one semi-transparent cylinder that are outside each other, or one 
is inside the other. 



IX. CASIMIR EFFECT OF A CYLINDER PARALLEL TO A PLANE 



In this section, we consider the Casimir interaction between a cylinder and a plane that are parallel to each 
other. Assume that the cylinder is represented by x 1 + y 2 = R 2 , —H/2 < z < H/2, and the plane is represented 
by x = L,—H/2 < y,z < H/2. The center of the cylinder is at O = (0,0,0) and the center of the plane is at 
O' = (L,0,0). 

In this case, we have to choose a different basis for plane waves. Choosing as the distinct direction, plane wave 
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basis are parametrized by (k y , k z ) £ M 2 , with 

reg 



tp° ut k ( x , fc) =e Ti V fe J- ~ k v x + ik yV+ ik z 



TE reg I reg 

A k fc ° Ut (x,fc) = . V X (f%* (x,fc) £ 



^fi^/k 2 . — k 2 x-\-ikuy-hik z z { • j „ -7 „ \ 

f v -l y " aT * (ik z e y — ik y e z ) , 



1.2 1 1,2 

TM reg 1 reg 

A fcHfc ; out (x,fc) = , V x V x ^(x.fcje, 



^ =e =FV*I-*S*-H*,»+*.« ((fc 2 + ^)e x ± kyy/kl - k*e y ± k zy /kl-k*e z ) 



k \ I ky -\- k 2 



Here k± = ^Jk 2 — fc 2 . 

Since the T-matrix for a cylinder and the T-matrix for a plane have been considered respectively in Section fVl and 
Section Hill let us now consider the translation matrices V and W. For our application, we only consider a translation 
by a vector L that is perpendicular to the z-axis. It is obvious that the translation matrices are diagonal in k z . In 
the scalar case, they are defined by 



8 fe >-L,fc) = £ W-L)^(x,fc), ( 86 ) 



Using (|78|) . we find that 



where 



Hence, 



r°° fik 

J —00 



(87) 



= f_^\n e —in<jib e —iy/k±—k2x+ikyy+ik z z 



COS0fe = - 



k 2 - kl 



7 2 + fc 2 , + fej, 



where 7 = ^/k 2 + fc 2 ■ 

For W, compare (1731 to ([57)1 , we find immediately that 



in<j> k J.\Jk\-k\L 



Hjkl-kl 



where 



cos <p k 



k 2 , - k 2 



k± 



sm <p k 
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Hence, 

W ky<n (Le x ,iZ) = —£= [ f^ +K ] e VW. 



Hence, for a scalar interaction between a cylinder and a plane, the Casimir interaction energy is given by 

TT POO pOO 11 

£cas = ^- / dU ^Tr In (1 - M(»0) , (88) 



where the trace Tr is given by (1801) . We have used the fact that the T, T and the translation matrices are diagonal in 
k z to take out the trace over k z . The matrix M is given by 



M n n' M = T nk * / V e in $ k T kvk * — c ™'<t>k c ?iJk 2 ,-klL 

'" 2n ' " Ik 2 , - kl 



V 

where T nkz is a component of the T-matrix of the cylinder derived in Section [Vj and T kykz is a component of the 
T-matrix of the plane derived in Section Hill with k± replaced by \Jky + k z . When the plane is subjected to Dirichlet 

or Neumann boundary conditions, T kykz = 1 or —1. Then 

/OO TT 

= ±T nk *^ n ,^(2Le x ). 

In imaginary frequency, 

M n , n ,(Le x ,i£) = ±T nk *K n+n ,{2 1 L). 
If general Robin boundary condition is imposed on the plane or the plane is semitransparent, we have 

n+n 



2 



Let 



Then 



k y = 7sinh(9. 



OO 



For a Robin plane, 



fk y k; = w-7COsh( 



u + 7 cosh 



For a semitransparent plane, 

A 



A + 27 cosh 6» 



In the formulas above, 7 = y/ k 2 + fcj, and k is related to £ by ([26 
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Next we turn to the electromagnetic case. As in the scalar case, the translation matrices are diagonal in k z . They 
are defined by 



A™>;-(x _ L, k) = £ «£ TE (-L)A™ : -(x, k) + V™' TE (-L)A™< -(x, *)) , 

n— — oo 

oo 

A™ : «*(x - L, k) = £ «£™(-L)A™; -(x, fc) + vSf™(-L)A™. reg (x, fc)) , 



(90) 



* TE. out / / . t 7 \ rr Z" 00 ^» / W TE,TE/ T \ A TE. out / I ,\ . ,,,TM,TE /T n a TM, out / / , \\ 

A n fc ; (x+L,fc)=#J_ -7^r[W kyn (liJA^^ (x,fc) + W^„ (L)A fcy ^ (x,fc)J, 

* TM, out / / . T , -i tj /" / TjrTE.TM/T \ * TE, out/ / ,\ , TI/ TM,TM/ T x A TM, out/ / , \ 

Kk z (x+L,fc)=#J_ —(W kyn (l)A kykz (x,k) + W kyr ; (L)A kykz (x,fc) 



The relation ^ implies that for Z = V or W, 

^TE.TM _ ^TM.TE ^TM,TM _ ^TE,TE 

For V, using (|85|) , we find that 



T /TE,TE/ t \ 



Q-n-A^K-^), 



<-l) n+i e 



71+1 il%4>k 



1 




1,2 

K z 


1 






kl 


k z cos 




+ 


z 



= — Le x 



t /TE.TM / r ->\ 



/il'Cg 



(S--A™: r ° S ) (-Le x ) 



=(-l)" +1 Q-n • 1 e -iy^fcF*+*Wfc.' ((fc^ + fca^ + kyy/kl-kfa + kjk\ - fcje, 



r--Le T 



(91) 



=(-l)" +1 e 



n+1 —in<f>k 



= (ism4> k e x -icoafaey) ■ ((fc 2 . + k 2 z )e x + k y ^k\ - k 2 y e y + k z ^jk\~ /^e 2 ) e'V^ fc ? L 



Y y ^ z 
zkku 



kj_ \J k y -\r k 2 
In imaginary frequency, 



T /TE,TE/ T 



^l 2 +k 2 y ( sji*+k?y+ky 



k 2 v + kl 7 



/-TE,TM ^ r _ ..^s fcy * K 



v n,k. (-Le x ,iQ=- 



\jl 2 +kl+ky 
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For W, since 
ism(f) k e x - i cos 4>ke y 



iky l y K x K v 
fci k\ 



ik y k ( y k y + k z k yy k l k l k z ykj_ k 2 \ k z ^k 2 _ k 2 I ^ ^ 



t. , lh.2 1 iL2 I iL /u2 \ U2 U lyi 1 £.2 / L . /pi 1.2 1 /jU2 1 jU2 /iU2 1 L2 

ft/J_ * / /i-y T rv^ \ \/ y «8 \/ y -2 / -L \ / y I \ \/ y ^ \/ y ^ 2 

Compare fl83| to (|5T]). we find that 

HJ-- + k 2 y Jk 2 + k 2 1 



Hence, for an electromagnetic Casimir interaction between a cylinder and a plane, the Casimir interaction energy is 
given by (|88|) . where 



In n , ( - = _ IjnK r dk v 1 ( k *\/^ + k2 y ~ ™ m ky \ ^ (k zy Jjl + kl ■ ;/■„,/, „ 



2 



^7m + fc y + V ~ iK rnk y k z ^^+k%J ^ -iK rn k y k z ^J"{^ + fc^ 

\/ "Trn 



n+n 



Here K m = n m {i^)^/c, 7 m = y + fc 2 and n m is the refractive index of the medium between the cylinder and the 
plane. Making a change of variables k y = 7 m sinh 6* , we find that 

M ,(if)=--T nk * [°° d9 - ( kzCOsh ® -iK m sinh9\ ^ kyk]t ( k z cosh9 -£/t TO sinh0\ 

2 J_ oc ^ sinh 2 (9 + fc 2 cosh 2 (9 \-iK m smh9 k z cosh9 J \-in m smh.9 k z cosh9 J 

v Jn+n')8 -2~j m L cosh 

Ac. c- 



If the plane is perfectly conducting, 
In this case, 

l„nnfc„ A 



1 

-1 



3G 



M n> „>(i£) = - ^T"* 5 * f J 1 ) / d0cosh([n + n']6V 



In general, T kj - is a diagonal matrix. If x > L is dielectric with permittivity e e and permeability /j, e , then 





A 2 


+ 7m 


sinh 2 £ 


1 - Me 7m coshfl 




Ae 2 


+ 7m 


sinh 2 £ 


1 + Me 7m COSh 


e ro \ 


111 


+ T 2 
1 ' m 


sinh 2 £» 


- e e 7 m cosh 9 



(92) 

p„. 1 / ^i 2 4- sinh 2 9 — cosh 



£ m \/ 7e + 7m sinh ^ + e e7m COsh 6» 



3G 



where j e = y/ k% + fcf, K e = n e (i£)£/c. 

The exact Casimir interaction energy for the cylinder-plane configuration has been considered in a number of 
works. For scalar interaction with Dirichlet or Neumann boundary conditions, or electromagnetic interaction with 
perfectly conducting boundary conditions, it has been considered in [3(| HH, Hy, H3, l4l| - |43l l7lj . The case where a 
magnetodielectric cylinder is in front of a magnetodielectric half-space was briefly discussed in [36| . In (38| , the Casimir 
interaction between a cylindrical plasma sheet and a dielectric half-space or a planar plasma sheet were treated as 
a scalar problem. It is interesting to note that the formula (jM)) for the scalar interaction between a cylinder and a 
plane is very similar to the formula obtained in (38j . 

In this section, we have also obtained the exact formula for the Casimir interaction energy for: 

• A Robin cylinder in front of a Robin plane. 

• A semitransparent cylinder in front of a semitransparent plane. 

• A Robin cylinder in front of a semitransparent plane. 

• A semitransparent cylinder in front of a Robin plane. 

X. CONCLUSIONS 

We have reported on the derivation of the Casimir interaction energy between two objects from the point of view 
of mode summation approach, for both scalar interactions and electromagnetic interactions. We consider the case 
where the two objects are exterior to each other, and the case where one object is inside the other. A closer scrutiny 
reveals that this approach actually has some flavor of the multiple scattering approach of [HI, [H, [3(| . The advantage 
over the multiple scattering approach is that we do not have to rely on the path integral quantization nor the Green's 
functions. In some sense, this provides a simpler prescription to obtain the Casimir interaction energy, and it is easier 
to generalize to other fields such as spinor fields and massive vector fields. 

In practice, all the difficulties in writing down the Casimir interaction energy lie in the computation of the T, T and 
the translation matrices, and the latter is usually more difficult. In the mode summation approach, one need to have 
a coordinate system for each of the objects where the wave equations are separable. In this paper, we consider planes, 
spheres and cylinders where the corresponding coordinate systems are readily available. After choosing a convenient 
basis, it is straightforward to compute the T and T matrices of an object by matching the boundary conditions. To 
compute the translation matrices, additional tools are needed, for which we have chosen the operator approach of 
48]. 

We discuss some situations that are of interest, such as a scalar field with Dirichlet, Neumann or Robin boundary 
conditions, or with a Dirac delta potential supported on the boundary of an object, and an electromagnetic field 
propagating in the presence of magnetodielectric object. This is not meant to be exhaustive but they are boundary 
conditions that have been explored so far. It is obvious that the mode summation approach can be straightforwardly 
generalized to other boundary conditions, which are only going to affect the T and T matrices of an object. 

In this paper, we have dealt exclusively with the zero temperature Casimir energy. Nevertheless, as mentioned at 
the end of Section [TH the extension to finite temperature is straightforward by using the Matsubara formalism. We 
only have to replace the formula (1151) for zero temperature Casimir energy by the formula (|23[) for finite temperature 
Casimir free energy. 

The results reported here do not only contain the known results about the sphere-sphere, cylinder-cylinder, sphere- 
plane and cylinder-plane configurations, but actually more. In particular, we have included results on scalar interaction 
between a Dirichlet/Neumann/Robin object and a semi-transparent object. The new results have been listed in each 
corresponding section. 
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Appendix A: Casimir interaction between multiple objects 

In this section, we explain how the mode summation approach in Section [TT] can be generalized to N objects, where 
N > 3. Let the N objects be labeled as Oi, O2, • ■ • , O^r. As in Section [Til we consider the following two cases 
separately: 

Case A The N objects are outside each other. 

Case B The objects Oi, O2, . . . , Oat_i are outside each other, but they are all inside the object On- 

Besides these two scenarios, there are other scenarios, such as object Oi is inside object O2, object O2 is inside 
object O3, and so on. The treatment of these scenarios are conceptually the same but technically more complicated. 
Therefore we only restrict ourselves here to the Case A and Case B mentioned above. For simplicity, we discuss only 
scalar fields. 

Let Oi be the coordinate origin of an appropriate coordinate system for object Oi, and let Xj be the position vector 
of a point with respect to Oi. If Lij is the vector from Oi to Oj, then Xi = Xj + Ly. Obviously, Ljj = — Ly. All the 
other notations are the same as in Section HU 



1. Case A 



The scalar field is represented by 



tp : 



/OO 
-00 



inside object Oi. In the region outside all the objects and close to object Oi 

Lp = 



/CO 
^^(a«V r a 7 (x . )W) + 6?V out (x . )a;)) f 
-OO 



(Al) 



(A2) 



As in Section HH the boundary conditions on the boundary of Oi give rise to two equations between A" 1 , a™ i and bf* 
Eliminating A?* gives a relation between af* and b" 1 : 



J i u i ' 



(A3) 



which defines T?*. On the other hand, the outgoing waves c/?° ut outside object Oi will propagate to other objects. In 
a vicinity of object Oj, where j ^ i, it contributes to tp™? via the translation relation: 

<f (x, , W ) = U aj , ai (Ly )^( Xj , U ) . (A4) 
In fact, (p r ° s is a superposition of the waves propagated from the outgoing waves close to other objects. Hence, 

(A5) 



N 



j = l Cti 

We are using the convention that U aitai (Tja) — 0. Let a and b be respectively the column matrices 



a 2 



b = 



b 2 



\h N J 



where ai is the infinite column vector with components a" 1 , and bi is the infinite column vector with components b° 
Let T and U be respectively the matrices given by 



T = 



/Ti ... \ 
T 2 ... 



\0 



T 



u = 



/ U12 

U21 



nJ 



Vin\ 
V2N 



\U N1 Um ... J 
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where Tj is an infinite diagonal matrix with (aj, aj)-element given by T- a< , and U,j is an infinite matrix with (cti,aj) 
given by U aiiCtj (Lji). (|A3I) and (IA5|) imply that 

b = -Ta = -TUb. 

As in Section HH this implies that the Casimir interaction energy is given by (|15l) . with 

M = -TU. 

In particular, when N = 2, 

-TiUu 
-T 2 U 2 i 



Hence, 



<!<•(.(] - K) = ( _ T l 2V2i T f 12 ) = <K (I - Tilic.T.U.,, ) . 



which recovers the result of Section HH 



2. Case B 



Inside the object 0^, I < i < JV — 1, the scalar field is still represented by (|AI j) . Outside the object Oat, the scalar 
field is represented as 



<P 



(A6) 



In the region outside the objects O^, I < i < N — 1 and inside O^r, the scalar field is represented as (|A2I) in a vicinity 
of Oj, 1 < j < N. The matrix Tj for the object Oi, 1 < i < N — I is defined by (|A3|) . For the object Otv, the matrix 
Tjv is defined by 



N N ' 



(A7) 



which is obtained by eliminating B N N from the boundary conditions on On- The translation formula (|A4j) still holds 
for 1 < i,j < iV — 1. Moreover, we have 



(A8) 



(A9) 



for 1 < i < JV - 1. These imply that for 1 < i < N - 1, 



3=1 Q i 



and 



AT-1 



i— 1 ttj 



Let 



a 2 



ajv_i 
V b N J 



b = 



( bl \ 
b 2 



bjv-i 



(A10) 



(All) 
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/Ti 
T 2 



T = 



\ 




Tjv/ 



/ 

Uai 



U = 









V o o 

OH), JM}, C^10|) and (jATTI) imply that 

b = -fa = -TUb 
Therefore, the Casimir interaction energy is given by (1151) . with 

M = -TO. 

When N = 2, 

I -T1V12 



V Wjvi W N2 



det (I - M) = det 
which recovers the result of Section ITT1 



-T 2 W. 



2 "21 



Ui,jv_i Yin \ 

U2.JV-I V2AT 



Vjv-i,jv 
Wvjv-i / 



det (I - T1V12T2W21) . 



Appendix B: Verification of identities 



First we verify the identity (|6ip . By definition, 

1 



X-i_ m (9 k , <f>k)-~X-i> m '(6k, 4>k) 



Let 



Then 



y/l(l + l)l'(l' + l) \sm 2 9 k 



-Yl-m(0k,4 > k)Yl> m '{0 k ,(i)k) + 



X;,- W (#fc, cf) k ) ■ X;> w /(gfe, fa) = ^ ^l",m"Yi" m "(6k,(f>k)- 



l"=0 



m" — — l 



de k 



(Bl) 



, -j> k 1 d9 k sin 

y/l(l + 1)1' (V + 1) Jo it) 



sin 2 



il,-fn(^*,^k)ll'm'(^,0fc) + 



dYi- m (6 k ,<pk) dY Vm >(0 k ,4>k) 



39 k 



d8 k 



Yi?, m „(6 k ,<t>k). 



For the integration with respect to 9 kl integration by parts gives: 

. Q dYi t - m (0 k ,(f>k) dYi> m >(9 k ,<f> k ) , 
d9 k sm 9 k — — Y v , m ,, {9 k ,<f> k ) 



d9 k 



sm V k Ou k 



d9 k 



Mi ^^-^tt4t (sme k dYl '- m a { ° k,(l>k) ) Y Vm ,{0k^k)Y{!, mll {9 k ^ k ) 



" sin 9 k d9 k 



d9 k sm9 k Y, 



-m\Vk,<Pk 



d9 k 

dY Vm ,{6 k ,4> k ) , &Y lt - m (9 k ,<l> k ) (a dYfi^ 

' * I' - 



d9 k 



d8 k 



'm'\Pk,9k 



08 k 



d9 k shx9 k I'Q' + 1) 



+ d9 k sin9 k (l(l + l) 



, Yi m (0 k ,(l) k )Y^ m :(9k,(t> k )Y l * m „(9 k ,^ k ) 

sin 9 k , 



Y l< - m (9 k ,<l> k )Y l , m ,(9 k ,c( >k )Y l % m „(9 k ,cl> k ) 



sin 2 9 k 



1 r 

2 Jo 



d9 k sm9 k [l"(l" + l)-—^- 
sm k 



Yi- m {9 k ,(j> k )Y Vml {9 k , <j) k )Y^, ml ,{6 k , 4> k ). 
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Using the fact that 

/■27T 

/ d<j) k Y L _ m (9 k , <j> k )Y Vm ,(9 k , <Ak)ij* m //(0fc) 4>k) 
Jo 

is nonzero only if ml = to + to", we find that 

id + hi + V(V + hi - + 1) r 2 * r 

Si»m// = — , = == == ' / d4> k / ^ fe sin^_ m (0 fe ,</» fe )y ;w (^,0 fe )^* m ,,(^,^) 
2 v //(/ + l)/'(Z' + l) Jo Jo 



+ !)/'(/' + 1) 



-(-l) r 



Z"(Z" + 1)- 1(1 + 1) -l'(l' + l) t « +1 /(2i + l)(2f' + l)(2i" + l) /'i /' 



4tt 



TO —TO TO 



This proves (joTj) . 

Next, we verify (|62l) . Since 



1 




hl)l'(l'- 


-1) 




ik 2 




y/l(l- 


-l)l'(V- 


-1) 




lk 2 




y/l(l- 


hl)l'{V- 


-1) 



m , ^ dY Vm ,(0 k ,4> k ) to ^ dYi_ m (9 k ,M 

m»t c% sm6» fc dd k 



VY Vm i{9k,<l>k) x VFi_ m (6> fc ,<fe)) • e fe 

V X {Yi, m ,(9 k ,(f> k )VY l ^ m {9 k ,4> k ) 



we find that 



27T />7T 



/, / d6 k sm9 k Xi- m (9 k , fa) ■ ( -r- x ~X-i> m >(9k,4>k) ) e l 



ii Jo 

d0fc / d^fe sin 6^, «L x 



ik 

ik 2 



ik-L 



JO 
k 



Y 



y/l(l + l)l'(l' + l) Jo 

Using the fact that 
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